Identification of bona fide RNA G-quadruplex binding proteins.

RNAs often accomplish their diverse functions through direct interactions with RNA-binding proteins (RBPs) in a sequence- and/or structure-dependent manner. RNA G-quadruplexes (rG4s) are unique secondary structures formed by guanine-rich RNA sequences which impact RNA function independently and in combination with RBPs. Efforts from several labs have identified dozens of rG4 specific RBPs (rG4BPs), although the research is still in the growing phase. Here we present methods for the systematic identification of rG4BPs using a pull-down approach that takes advantage of the chemical modification of guanine bases. This allows abolishing the rG4 structures while still maintaining the base composition intact, and hence helps in recognizing true rG4BPS (in contrast to G-rich motif binders). In combination with other biochemical assays, such an approach can be efficiently used for the identification and characterization of bona fide rG4BPs.

Structure of native four-repeat satellite III sequence with non-canonical base interactions.

Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5'-ATGGA ATGGA ATGGA ATGGA-3' using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05-0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G-G intercalation sandwiched by sheared A-G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A-T/T-A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.

Formation and Nanomechanical Properties of Silver-Mediated Guanine DNA Duplexes in Aqueous Solution.

Silver cations can mediate base pairing of guanine (G) DNA oligomers, yielding linear parallel G-Ag+-G duplexes with enhanced stabilities compared to those of canonical DNA duplexes. To enable their use in programmable DNA nanotechnologies, it is critical to understand solution-state formation and the nanomechanical stiffness of G-Ag+-G duplexes. Using temperature-controlled circular dichroism (CD) spectroscopy, we find that heating mixtures of G oligomers and silver salt above 50 °C fully destabilizes G-quadruplex structures and converts oligomers to G-Ag+-G duplexes. Electrospray ionization mass spectrometry supports that G-Ag+-G duplexes form at stoichiometries of 1 Ag+ per base pair, and CD spectroscopy suggests that as the Ag+/base stoichiometry increases further, G-Ag+-G duplexes undergo additional morphological changes. Using liquid-phase atomic force microscopy, we find that this excess Ag+ enables assembly of long fiberlike structures with ∼2.5 nm heights equivalent to a single DNA duplex but with lengths that far exceed a single duplex. Finally, using the conditions established to form single G-Ag+-G duplexes, we use a surface forces apparatus (SFA) to compare the solution-phase stiffness of single G-Ag+-G duplexes with dG-dC Watson-Crick-Franklin duplexes. SFA shows that G-Ag+-G duplexes are 1.3 times stiffer than dG-dC duplexes, confirming gas-phase ion mobility spectrometry measurements and computational predictions. These findings may guide the development of structural DNA nanotechnologies that rely on silver-mediated base pairing.

Structure of a DNA G-quadruplex that Modulates SP1 Binding Sites Architecture in HIV-1 Promoter.

Nucleic acid sequences containing guanine tracts are able to form non-canonical DNA or RNA structures known as G-quadruplexes (or G4s). These structures, based on the stacking of G-tetrads, are involved in various biological processes such as gene expression regulation. Here, we investigated a G4 forming sequence, HIVpro2, derived from the HIV-1 promoter. This motif is located 60 nucleotides upstream of the proviral Transcription Starting Site (TSS) and overlaps with two SP1 transcription factor binding sites. Using NMR spectroscopy, we determined that HIVpro2 forms a hybrid type G4 structure with a core that is interrupted by a single nucleotide bulge. An additional reverse-Hoogsteen AT base pair is stacked on top of the tetrad. SP1 transcription factor is known to regulate transcription activity of many genes through the recognition of Guanine-rich duplex motifs. Here, the formation of HIVpro2 G4 may modulate SP1 binding sites architecture by competing with the formation of the canonical duplex structure. Such DNA structural switch potentially participates to the regulation of viral transcription and may also interfere with HIV-1 reactivation or viral latency.

Best method to determine DNA G-quadruplex folding: The 1H-13C HSQC NMR experiment.

NMR spectroscopy is the major method for G-quadruplex structure determination under physiologically relevant solution conditions. Unlike duplex B-DNA, in which all nucleotides adopt an anti glycosidic conformation, the core tetrad-guanines in a G-quadruplex can adopt anti or syn glycosidic conformation depending on the folding structure. An experimental method that can clearly and unambiguously determine syn and anti tetrad-Gs in a G-quadruplex is highly desirable and necessary. In the present study, we exploit the advantages of the 1H-13C HSQC experiment to determine tetrad-G's glycosidic conformation and thus folding topology of G-quadruplexes. We use several examples to demonstrate the clear and straightforward determination of the guanine glycosidic conformations and G-quadruplex folding structures. Moreover, 1H-13C HSQC data can readily identify adenine H2 resonances as well as determine unusual syn conformation in loop and flanking sequences, a challenging task by standard 2D NOESY.

Base flipping mechanism and binding strength of methyl-damaged DNA during the interaction with AGT.

Methyl damage to DNA bases is common in the cell nucleus. O6-alkylguanine-DNA alkyl transferase (AGT) may be a promising candidate for direct damage reversal in methylated DNA (mDNA) at the O6 point of the guanine. Indeed, atomic-level investigations in the contact region of AGT-DNA complex can provide an in-depth understanding of their binding mechanism, allowing to evaluate the silico-drug nature of AGT and its utility in removing methyl damage in DNA. In this study, molecular dynamics (MD) simulation was utilized to examine the flipping of methylated nucleotide, the binding mechanism between mDNA and AGT, and the comparison of binding strength prior and post methyl transfer to AGT. The study reveals that methylation at the O6 atom of guanine weakens the hydrogen bond (H-bond) between guanine and cytosine, permitting for the flipping of such nucleotide. The formation of a H-bond between the base pair of methylated nucleotide (i.e., cytosine) and the intercalated arginine of AGT also forces the nucleotide to rotate. Following that, electrostatics and van der Waals contacts as well as hydrogen bonding contribute to form the complex of DNA and protein. The stronger binding of AGT with DNA before methyl transfer creates the suitable condition to transfer methyl adduct from DNA to AGT.

Quantitative Structure-Activity Relationship in the Series of 5-Ethyluridine, N2-Guanine, and 6-Oxopurine Derivatives with Pronounced Anti-Herpetic Activity.

A quantitative analysis of the relationship between the structure and inhibitory activity against the herpes simplex virus thymidine kinase (HSV-TK) was performed for the series of 5-ethyluridine, N2-guanine, and 6-oxopurines derivatives with pronounced anti-herpetic activity (IC50 = 0.09 ÷ 160,000 µmol/L) using the GUSAR 2019 software. On the basis of the MNA and QNA descriptors and whole-molecule descriptors using the self-consistent regression, 12 statistically significant consensus models for predicting numerical pIC50 values were constructed. These models demonstrated high predictive accuracy for the training and test sets. Molecular fragments of HSV-1 and HSV-2 TK inhibitors that enhance or diminish the anti-herpetic activity are considered. Virtual screening of the ChEMBL database using the developed QSAR models revealed 42 new effective HSV-1 and HSV-2 TK inhibitors. These compounds are promising for further research. The obtained data open up new opportunities for developing novel effective inhibitors of TK.

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA.

Proton transfer across hydrogen bonds in DNA can produce non-canonical nucleobase dimers and is a possible source of single-point mutations when these forms mismatch under replication. Previous computational studies have revealed this process to be energetically feasible for the guanine-cytosine (GC) base pair, but the tautomeric product (G[Formula: see text]C[Formula: see text]) is short-lived. In this work we reveal, for the first time, the direct effect of the replisome enzymes on proton transfer, rectifying the shortcomings of existing models. Multi-scale quantum mechanical/molecular dynamics (QM/MM) simulations reveal the effect of the bacterial PcrA Helicase on the double proton transfer in the GC base pair. It is shown that the local protein environment drastically increases the activation and reaction energies for the double proton transfer, modifying the tautomeric equilibrium. We propose a regime in which the proton transfer is dominated by tunnelling, taking place instantaneously and without atomic rearrangement of the local environment. In this paradigm, we can reconcile the metastable nature of the tautomer and show that ensemble averaging methods obscure detail in the reaction profile. Our results highlight the importance of explicit environmental models and suggest that asparagine N624 serves a secondary function of reducing spontaneous mutations in PcrA Helicase.

Strengthened cooperativity of DNA-based cyclic hydrogen-bonded rosettes by subtle functionalization.

Cooperative effects cause extra stabilization of hydrogen-bonded supramolecular systems. In this work we have designed hydrogen-bonded rosettes derived from a guanine-cytosine Janus-type motif with the aim of finding a monomer that enhances the synergy of supramolecular systems. For this, relativistic dispersion-corrected density functional theory computations have been performed. Our proposal involves a monomer with three hydrogen-bonds pointing in the same direction, which translates into shorter bonds, stronger donor-acceptor interactions, and more attractive electrostatic interactions, thus giving rise to rosettes with strengthened cooperativity. This newly designed rosette has triple the cooperativity found for the naturally occurring guanine quadruplex.

Effect of Guanine Adduct Size, Shape, and Linker Type on the Conformation of Adducted DNA: A DFT and Molecular Dynamics Study.

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 µs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), ß (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti-to-syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with ß-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

Investigating the Spectroscopy of the Gas Phase Guanine-Cytosine Pair: Keto versus Enol Configurations.

We report on a vibrational study of the guanine-cytosine dimer tautomers using state-of-the-art quasiclassical trajectory and semiclassical vibrational spectroscopy. The latter includes possible quantum mechanical effects. Through an accurate comparison to the experimental spectra, we are able to shine a light on the hydrogen bond network of one of the main subunits of DNA and put the experimental assignment on a solid footing. Our calculations corroborate the experimental conclusion that the global minimum Watson-and-Crick structure is not detected in the spectra, and there is no evidence of tunnel-effect-based double proton hopping. Our accurate assignment of the spectral features may also serve as a basis for the development of precise force fields to study the guanine-cytosine dimer.

8-oxoguanine riboswitches in bacteria detect and respond to oxidative DNA damage.

Riboswitches rely on structured aptamer domains to selectively sense their target ligands and regulate gene expression. However, some riboswitch aptamers in bacteria carry mutations in their otherwise strictly conserved binding pockets that change ligand specificities. The aptamer domain of a riboswitch class originally found to selectively sense guanine forms a three-stem junction that has since been observed to exploit numerous alterations in its ligand-binding pocket. These rare variants have modified their ligand specificities to sense other purines or purine derivatives, including adenine, 2'-deoxyguanosine (three classes), and xanthine. Herein, we report the characteristics of a rare variant that is narrowly distributed in the Paenibacillaceae family of bacteria. Known representatives are always associated with genes encoding 8-oxoguanine deaminase. As predicted from this gene association, these variant riboswitches tightly bind 8-oxoguanine (8-oxoG), strongly discriminate against other purine derivatives, and function as genetic "ON" switches. Following exposure of cells to certain oxidative stresses, a representative 8-oxoG riboswitch activates gene expression, likely caused by the accumulation of 8-oxoG due to oxidative damage to G nucleobases in DNA, RNA, and the nucleotide pool. Furthermore, an engineered version of the variant aptamer was prepared that exhibits specificity for 8-oxoadenine, further demonstrating that RNA aptamers can acquire mutations that expand their ability to detect and respond to oxidative damage.

Reaction of human telomeric unit TTAGGG and a photoactivatable Pt(IV) anticancer prodrug.

The interaction of a photoactivatable diazidodihydroxido Pt(IV) prodrug, trans,trans,trans-[Pt(N3)2(OH)2(py)2] (py = pyridine; 1), with a hexamer straight human telomeric DNA unit sequence (5'-T1T2A3G4G5G6-3', I) upon light irradiation was investigated by electrospray ionization mass spectroscopy (ESI-MS). In the primary mass spectrum, two major mono-platinated I adducts with the bound Pt moieties, trans-[PtII(N3)(py)2]+ (1') and trans-[PtII(py)2]2+ (1''), respectively, were detected. It is rare to observe such high abundance and nearly equal intensity platinated DNA adducts formed by these two PtII species because 1' is usually the only major reduced Pt(II) species produced by the photodecomposition of complex 1 in the presence of DNA while 1'' was rarely detected as the major reduced PtII species reported previously. Subsequent tandem mass spectrometric analysis by collision-induced dissociation (CID) showed that in the former adduct {I + 1'}2+, G6 and A3 were the platination sites. While in the latter adduct {I + 1''}2+, a potential intrastrand crosslink was speculated after G4 and G6 sites were identified. Additionally, other minor platinated adducts like di-platinated I adduct by 1' with platination sites at G4 and G6 and mono-platinated I adducts containing base oxidation were also detected by mass spectrometry. Due to the rich guanines and their sensitivity to oxidation, the oxidation induced by 1 most probably occurred at guanine. The oxidation adducts were proposed as 8-hydroxyl guanine, spiroiminodihydantoin (Sp), 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapyG), 5-guanidinohydantoin (Gh), and/or dehydroguanidinohydantoin (DGh) referring to previous reports. The obtained results provide useful chemical information about the photoreaction between photoactivatable Pt(IV) anticancer prodrugs and human telomeric DNA. Such special damages of Pt(IV) prodrugs on human telomeric DNA implicate its active role in the mechanism of Pt(IV) prodrugs and further support the unique sequence-dependent photointeraction profile of complex 1 reacting with DNA.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA.

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations.

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

Steric hindrance and structural flexibility shape the functional properties of a guanine-rich oligonucleotide.

Ligand/protein molecular recognition involves a dynamic process, whereby both partners require a degree of structural plasticity to regulate the binding/unbinding event. Here, we present the characterization of the interaction between a highly dynamic G-rich oligonucleotide, M08s-1, and its target protein, human α-thrombin. M08s-1 is the most active anticoagulant aptamer selected thus far. Circular dichroism and gel electrophoresis analyses indicate that both intramolecular and intermolecular G-quadruplex structures are populated in solution. The presence of thrombin stabilises the antiparallel intramolecular chair-like G-quadruplex conformation, that provides by far the main contribution to the biological activity of the aptamer. The crystal structure of the thrombin-oligonucleotide complex reveals that M08s-1 adopts a kinked structural organization formed by a G-quadruplex domain and a long duplex module, linked by a stretch of five purine bases. The quadruplex motif hooks the exosite I region of thrombin and the duplex region is folded towards the surface of the protein. This structural feature, which has never been observed in other anti-exosite I aptamers with a shorter duplex motif, hinders the approach of a protein substrate to the active site region and may well explain the significant increase in the anticoagulant activity of M08s-1 compared to the other anti-exosite I aptamers.

Development of Mn2+-Specific Biosensor Using G-Quadruplex-Based DNA.

Metal ions are used in various situations in living organisms and as a part of functional materials. Since the excessive intake of metal ions can cause health hazards and environmental pollution, the development of new molecules that can monitor metal ion concentrations with high sensitivity and selectivity is strongly desired. DNA can form various structures, and these structures and their properties have been used in a wide range of fields, including materials, sensors, and drugs. Guanine-rich sequences respond to metal ions and form G-quadruplex structures and G-wires, which are the self-assembling macromolecules of G-quadruplex structures. Therefore, guanine-rich DNA can be applied to a metal ion-detection sensor and functional materials. In this study, the IRDAptamer library originally designed based on G-quadruplex structures was used to screen for Mn2+, which is known to induce neurodegenerative diseases. Circular dichroism and fluorescence analysis using Thioflavin T showed that the identified IRDAptamer sequence designated MnG4C1 forms a non-canonical G-quadruplex structure in response to low concentrations of Mn2+. A serum resistance and thermostability analysis revealed that MnG4C1 acquired stability in a Mn2+-dependent manner. A Förster resonance energy transfer (FRET) system using fluorescent molecules attached to the termini of MnG4C1 showed that FRET was effectively induced based on Mn2+-dependent conformational changes, and the limit of detection (LOD) was 0.76 µM for Mn2+. These results suggested that MnG4C1 can be used as a novel DNA-based Mn2+-detecting molecule.

RNA structural probing of guanine and uracil nucleotides in yeast.

RNA structure can be essential for its cellular function. Therefore, methods to investigate the structure of RNA in vivo are of great importance for understanding the role of cellular RNAs. RNA structure probing is an indirect method to asess the three-dimensional structure of RNA by analyzing the reactivity of different nucleotides to chemical modifications. Dimethyl sulfate (DMS) is a well-established compound that reports on base pairing context of adenine (A) and cytidine (C) in-vitro and in-vivo, but is not reactive to guanine (G) or uracil (U). Recently, new compounds were used to modify Gs and Us in plant, bacteria, and human cells. To complement the scope of RNA structural probing by chemical modifications in the model organism yeast, we analyze the effectiveness of guanine modification by the glyoxal family in Saccharomyces cerevisiae and Candida albicans. We show that within glyoxal family of compounds, phenylglyoxal (PGO) is the best guanine probe for structural probing in S. cerevisiae and C. albicans. Further, we show that PGO treatment does not affect the processing of different RNA species in the cell and is not toxic for the cells under the conditions we have established for RNA structural probing. We also explore the effectiveness of uracil modification by Cyclohexyl-3-(2-Morpholinoethyl) Carbodiimide metho-p-Toluenesulfonate (CMCT) in vivo and demonstrate that uracils can be modified by CMCT in S. cerevisiae in vivo. Our results provide the conditions for in vivo probing the reactivity of guanine and uracil nucleotides in RNA structures in yeast and offer a valuable tool for studying RNA structure and function in two widely used yeast model systems.

Theoretical insights into the reaction mechanism of hydroxyl radicals and guanine in G-quadruplex DNA.

A DFT investigation was performed to illuminate the obscure mechanism of hydroxyl radical (OH˙) and guanine in G-quadruplex by mapping the energy profiles for both addition and hydrogen abstraction reactions. Results revealed that in G-quadruplex, the electrophilic attack of OH˙ to C8 (G) leading to 8-oxoG is the most energetically favorable course, where direct hydrogen abstraction from N2 of G to furnish neutral radicals could compete with that. Although the addition of OH˙ to C4 and C5 positions could provide stable OH-adducts, the subsequent dehydration of C4-OH adduct and hydrogen transfer of C5-OH adduct, which is a prerequisite for neutral radical formation, is rate-limited due to the high barrier manifesting the inaccessibility for these courses. Intriguingly, the identity of the decisive neutral radical was confirmed to be G(N2-H)˙ rather than the familiar G(N1-H)˙, where the hydrogen bond plays significant roles by blocking tautomerizations.

[Oxidative Probing of the G4 DNA Structure by Znp1 Porphyrin within Sequences of MYC and TERT Promotors].

The formation of G4 structures in a DNA double helix competes with the complementary strand interaction. The local environment in DNA can change equilibrium of G4 structures, which are studied on single-stranded (ss) models by classical structural methods. A relevant task is to develop methods for detecting and localizing G4 structures in extended native double-stranded (ds) DNA in the promoter regions of the genome. The ZnP1 porphyrin derivative selectively binds to G4 structures and leads to photo-induced oxidation of guanine in ssDNA and dsDNA model systems. We have shown the oxidative effect of ZnP1 on native sequences of MYC and TERT oncogene promoters, which can form G4 structures. Single-strand breaks in the guanine-rich sequence because of ZnP1 oxidation and subsequent cleavage of the DNA strand with Fpg glycosylase have been identified and assigned to the nucleotide sequence. The detected break sites have been shown to correspond to sequences capable of forming G4 structures. Thus, we have demonstrated the possibility of using porphyrin ZnP1 for the identification and localization of G4 quadruplexes in extended regions of the genome. Here we have shown the novel data on a possibility of folding G4 structures in the presence of complementary strand in native DNA double helix.